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Abstract. The electronic entanglement between two free atoms initially at rest is obtained 
including the effects of photon recoil, for the case when quantum dispersion can be neglected 
during the atomic excited-state lifetime. Different from previous treatments using common or 
statistically independent reservoirs, a continuous transition between these limits is observed, 
that depends on the inter-atomic distance and degree of localization. The occurrence of 
entanglement sudden death and birth as predicted here deviates from the case where the inter- 
atomic distance is treated classically by a static value. Moreover, the creation of a dark state is 
predicted, which manifests itself by a stationary entanglement that even may be created from 
an initially separable state. 
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1. Introduction 

If a quantum system is brought into contact with another system with much larger mode 
space, the system suffers decoherence HI, which in Born-Markov approximation results as 
an exponential decay of its coherences. Within this paradigm one may conclude that also 
entanglement J2l should naturally suffer a degradation with a similar behavior: an asymptotic 
decay to zero. However, recently it has been shown that many exceptions from this expected 
behavior exist, that show a vanishing of entanglement at finite time, including its possible 
reappearance at later times. These features have been denoted as entanglement sudden death 
(ESD) and birth (ESB), respectively (10. 

The occurrence of ESD and/or ESB has created increasing interest in studying the 
dynamics of entanglement of various types of bipartite systems, where the concurrence ifTOl 
can be used as a measure of entanglement. These systems include among others, systems 
of continuous variables JED EL atoms in cavities fTTl IT3l - iT5ll . cavity fields lfT6l[P71l . spin 
chains |[T8l[T9ll . and atoms in free space ||20] - |231 . Also non-Markovian effects have been 
studied for the dynamics of bipartite entanglement Il24l - [2"6ll . Furthermore, experiments have 
been performed to model [|27ll or directly measure ll28l the dynamics of entanglement. 

To best of our knowledge, at present, theoretical analysis has been performed only within 
the framework of a master equation in the Born-Markov approximation, where the atoms are 
supposed to have classical and fixed positions 129113011 . In such a treatment 112111 the dynamics 
of the concurrence has been shown to reveal in general two different time scales, one given 
by the spontaneous decay rate and the other by a distance-dependent collective decay rate. 
Furthermore, the existence of a sequence of more than one ESD and one ESB has been stated. 
In the limit of infinite inter-atomic distance, the well-known results for two spin systems in 
statistically independent reservoirs result [7J. However, in such a treatment the influence of 
the relative quantized motion of the atoms is discarded. 

Physically, this may only be realized by situating the two atoms in independent traps and 
in the limit where the trap frequency is CO ^> E r /h, where E t is the recoil energy of the atoms. 
In this Lamb-Dicke limit the photon recoil is received by the trapping mechanism (Mossbauer 
effect) and the atomic motion is unaffected. In all other cases, the photon recoil affects the 
motion of the two atoms and in this way the dynamics of the electronic entanglement between 
the atoms may be modified. This is most relevant in the limit of free atoms. There the relative 
quantized motion has to be included, in order to allow the distinguishability between the atoms 
for properly addressing the issue of their mutual entanglement ll23l . We will briefly review 
the arguments given in Ref. Il23l : 

Given two identical two-level atoms, initially at rest and distinguishable by their differing 
positions, their distinguishability should be maintained at least during the atom's excited-state 
lifetime To = yo being the natural linewidth of the electronic transition. This condition 

requires that the quantum dispersion of the relative-position wave-packet be sufficiently weak, 
so that during To the rms spread be 

Ar < r, (1) 

where f is the mean inter-atomic distance. Due to quantum dispersion, the initial rms spread 
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Arj is enlarged during the interval Tq to 




where the dispersion length is / d = ^/h/jom with h being Planck's constant and m the atomic 



mass. 



Equation © together with condition © leads to the condition 



'»^»Ao^, (3) 
where Ao is the wavelength of the electronic transition and the recoil energy is E r = (M ) 2 /2m 




with &o = 2% /Xq. As typically E r <C hy^, the minimum distance and rms spread may still 
be much smaller than Xq. However, a zero rms spread is not permitted, as then quantum 
dispersion would rapidly render the atoms indistinguishable. Thus, to be consistent with the 
requirement of distinguishability during Xq, a finite initial spread Ar± within the limits © is 
required. 

This latter requirement implies the quantized treatment of the inter-atomic motion, which 
includes the recoil effects of photon emissions. In consequence, in this paper we present a 
fully quantized treatment of both the electronic and relative motion. By use of a Wigner- 
Weisskopf approach it is shown how the dynamics of the entanglement between the atoms is 
affected by the photon recoil during spontaneous emissions. 

The outline of the paper is as follows: In Sec. [2] the interaction between the two atoms 
and the em field is treated in the Wigner-Weisskopf approach and the dynamics of the relevant 
probability amplitudes is solved. Using these solutions, in Sec. Ill the electronic reduced 
density matrix is obtained and then used in Sec. IV to derive the concurrence as a measure of 
entanglement. Finally, in Sec. V a summary and conclusions are given. 

2. Dynamics of two two-level atoms interacting with the electromagnetic field 

2.1. Hamilton operator 

The system consists of two atoms of mass m and with inter-atomic distance vector r and 
corresponding relative momentum p. Each atom is composed of two electronic levels with 
transition frequency C0q, described by the atomic pseudo spin operator S a (a = ±1). Both 
atoms interact with the em field via their electronic transition dipole moment d. 
The complete Hamilton operator for this system is given by 




(4) 



Here d-^ a is the annihilation operator for a photon in the plane-wave mode defined by the 
wave vector k and the polarization a. The interaction between atoms and em field reads 




(5) 
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with a being the polarization unit vector and with the vacuum electric -field strength 



hck 

167T 3 £o' 



(6) 



(7) 



2.2. Equations of motion for the probability amplitudes 

Assuming an initial vacuum state for the em field and an initial electronic state of the atoms 
consisting of a superposition of only | W) and | 44) > the general form of the wave function for 
the system is given by 



\ x ¥(t)}= J d 3 r\r) rel ® ^\j/ g (r,t)\g} + ^(7,^) 

+ I ^LLv^fcOk.M) 

J a a 

+ [d 3 k [ d 3 k' £ \fy G £ a ,(r,t)\k,o-,k',o') 

J J n ' 



Here the electronic-photonic states are defined as 



Is) 

\e) 

\+;k,a) 
\-;k,a) 
\k, o;k', o') 



44) el® | Vac) em, 
tt)d®|vac)em, 

4t)el®4al VaC ^ e 
at 



t4)el®^ a |vac)em, 
44)el®4cr4a'l VaC ^ 



(8) 

(9) 
(10) 
(ID 

(12) 

(13) 



where the state (IT3T ) is not normalized to unity in order to contain the bosonic enhancement 
when two identical photons are emitted. 

With the general form © the Laplace-transformed Schrodinger equation, 

ihs\W(s)) = ih\<i>} +ft\¥(s)), (14) 

has to be solved with the initial condition |*P(f = 0)) = |4>). The initial state |4>) is analogously 
expanded as Eq. ([8]) with the initial probability amplitudes being <j> g (r), (j> e (r), whereas 
ty a .% a (f) = 0£ a .%r a i(r) = because of the initial absence of photons. The corresponding 
Laplace-transformed coupled equations for the probability amplitudes read then 

(is + (Oo)\j/ g (r,s)=i<j) g (f), (15) 

(w-fl*)V.ft') = '*e(7) + /^EE^a^fr'y^. (16) 



X 



-iak! -r/2 



(17) 
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{is+COo-ck-ck') n,a»,o^ S )=L ^a^^)^ 2 

a ' 

+ 4,a^,o^y^ /2 ] ■ (18) 
Here it was assumed that the resulting spectral width of the electronic resonances is much 
larger than the Doppler shifts induced by the slow relative motion. Under this assumption the 
kinetic energy of this relative motion can be safely neglected. 

Note, that with the appearance of the state |r) re i in Eq. ([8]) a full quantized treatment 
of the relative distance between the atoms is warranted and included by the use of the 
correspondingly defined amplitudes. As a consequence, the neglection of the kinetic energy 
in the equations of motion for the amplitudes, cf. Eqs. (fT5l) - (fT8l) . is equivalent to the Raman- 
Nath approximation in atomic -beam diffraction. There the kinetic energy is neglected, but the 
photon recoil still generates orders of diffraction of the passing atoms which correspond to 
dynamical changes of their atomic momenta. 

2.3. Solution of the probability amplitudes 

From Eq. (fT5l) is can be readily seen that the ground-state amplitude has the simple solution 

x i r g {r,t) = <t> g {r)e i <* t . (19) 

To obtain the solution for the excited-state amplitude \j/ e , we first insert Eq. (PT8l) into Eq. (fTTT) 
to obtain 



(is-ck) W a -X,a^ s ) = K L^^^ e 



-iak-7/2 



+ 1**1.7. 



is + 0)q — c(k + k') 

+1**1* 



k'-J 



2 KV/ _/ K"— 

k',o' la 



x 



(20) 



~f is+ CQo — c(k + k') 

W a . X , o/ (7,s)e'^-^ 
Assuming now the form 

and assuming that we may expand the required solution % .? in powers of the atom-photon 
interaction rate k? „, 

we obtain from the first-order terms of Eq. (1201) the condition 

(23) 



a;k,a is — ck 

and likewise from the third-order terms the condition 



\ K k'a' 



is + (Oq -c(k + k') 
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Inserting then Eq. (|2TI) into Eq. (fT6l) and using the expansion (|22l) . the excited state 
amplitude is obtained as 

xj/ e (r,s) =f(r,s)(j> e (r), (25) 

where the function is defined as 

with the complex- valued frequency shift being 

Q(r, s) = Q (1) (r, 5) + ft (3) (r, j) + . . . (27) 
The first-order term of this frequency shift results from and reads 

^ 1 )(r,,) = 2/j 3 fc£i^, (28) 



a is — ck 



whereas the third-order term results from %^ as 



^' ^ I ^ is — ck J ^ is + COq — ck — ck' 

" (29) 

We note that each occurrence of the factor | | ^ in Eqs (1281) and (l29l) will result as the 
natural line width 70 of the atomic electronic transition, which is defined as 

*=3^ <30) 

Therefore QW(r,j) °= 70 and Q.( 3 \r,s) °= 7q, where 70 <C (% so that the pole in Eq. (|26|) is 
located very near to 5 = — zft)0 and Q. represents a relative displacement of this pole of the 
order of e = 7o/tt>o <C 1 . Thus we may approximate this frequency shift by taking the limiting 
value 

Q(r)« lim S2(r,j) = lim QW(r,s) + (31) 



The limiting value is obtained by evaluating the residues of expression (1281) as 
lim/^^ QW(r, = — /70, so that we may use the approximation 

(32) 

In the time domain, the function f(r,t) can thus be approximated as 

f(r, t) « /(f) = exp [-(ifflo + 70)?] , (33) 

being independent of the inter-atomic distance vector r. Correspondingly the amplitude (|2TI) 
can be consistently approximated as 

■£* g—iak-r/2 

^, a (r,s)= l l_ ck fm(f). (34) 

Equations (fT9l) , (1231) and (1341) together with the approximation (|33l) constitute the required 
information to evaluate the reduced electronic density operator of the system. 
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3.1. Reduced electronic density matrix 

The reduced density operator of the electronic systems of the two atoms reads 
p(t)= fd 3 rTr em (r\xif(t))(x{f(t)\r). 



(35) 

In the standard basis {| ft), 1 | | H)} it can written as the reduced electronic density 
matrix 



P(0 







Pe,e 
C7+ 

o <j; 









Pe,g 

o+_ 

<7__ 
V Pe,g " P 8 ,8 J 

where the matrix elements are given by 
Pe,eit) = J d 3 r| !//,(?, f)| 2 , 

Pe, g (t) = J J 3 rv/ e (r,5)v/|(r,0, 

o a , b (t) = fd 3 r j d 3 kY,%.i a (7j)V* b . t(7 (lt) 
J J a 

and due to the unit trace, 

Pg,g( t ) = 1 -PeA t )- E °a,a{t)- 

a=± 



(36) 



(37) 
(38) 
(39) 

(40) 



For the evaluation of these matrix elements we will assume initial conditions where the 
relative motion of the atoms is factorized from the electronic state, i.e. 

|0 e (r)| 2 =p e , e (0)w(r), (41) 
\(j) g (r)\ 2 = Pg , g (0)w(r), (42) 

<j> e (^r g (r) = Pe, g ^Mr), (43) 
where w(r) is the probability density for the inter-atomic distance vector r. 

Given the approximation (|33l) . the solutions for the excited state population and of the 
outer off-diagonal matrix element become then 

p e , e (0=p e , e (0)e- 2 ^, (44) 
p e , s (t)=p e , 8 (0)e- 2 '^-™. (45) 

Due to the two possible decay channels | ff) | |!) and | |t) I It). tne decay of the 
excited state population occurs at twice the decay rate of a single atom. 

3.2. Density matrix elements o U) b 

From Eq. d39l) the Laplace transformed diagonal elements 0± ) ±(s) can be shown to be 
identical and of the form 



CT±,±(*) = p e , e (0) / ds'f(r : s-s')f*(r,s')I(s-s' : s') : 
Jo 



(46) 
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where the averaging is done over the probability density w(r) and the integral / is given by 

(s>p) 

*>p) = f d3 % { JwLi P y (47) 

After performing the sum over polarizations and integrating over spherical angles of the wave 
vector, this integral becomes^! 

I(s,p) = — J% / d(0 7 r . -. (48) 

IkG^JO (co-is)((0 + ip) 

Using Eqs (1331) and (1481) . the matrix elements (1461) can be rewritten as 

°±M S ) = ^r^Pe,(0) f°<W f ' ds'g(a),s-s')g*(co,s'), (49) 
2na>Q Jo Jo 

where the new function g is given by 

g ((0,s) = ^^ = - . . . - , -, (50) 

v ' co-is (s + ico)(s + icoo + Yo) 



where the definition (|26l) together with Eq. (1321) have been used. Thus, in the time domain the 
diagonal matrix elements are 

O±,±(0 = -^p e , e (0) / d(OCQ 3 \g(CQ,t)\ 2 . (51) 
ItzcOq Jo 

The Laplace transformed centered off-diagonal matrix elements can by obtained 
analogously by inserting Eq. (1341) into Eq. d39l) . Both elements result identical and read 



o±, T {s) = p e>e (0) jf ds'f(r,s-s')f*(r,s% (s -s',s',x r ), (52) 

where the retardation time of light propagation between the atoms is x r = r/c. The integral is 
defined as 

\ K ickz cos 6 

^■^ = / A | feg (sip) - <53) 

with £ = sin 2 Z(r, J), assuming the induced electrical dipole moments of the two atoms to 
be parallel. After summing over polarizations and integrating over the spherical angles of the 
wavevector, this integral becomes 

h(s,p,T) = — f \ dcoco 3 - \\ . . , (s>p), (54) 

where the dissipative part of the dipole-dipole interaction pattern is given by [|29l[30l 

(55) 



( ^ 3 



^ , cosx smx \ y. sin* 

(3^-2)(— --r)+^ 



% Note that the divergence of this integral is removed by the fact that it is integrated over s and p, as given in 
Eq. (®. 



Dynamics of entanglement between two free atoms with quantized motion 9 

Inserting Eq. (|54l ) into ( 1521 and using the approximation (1331 the off-diagonal elements 
result as 

0±, T (*) = ^^V(°)/ dCOC0 3 ll^((Or/c) 



X 



/ Jj'g(fi),s-s')g*(fi),/). (56) 
Jo 



Thus, in the time domain the matrix elements (l56l) become 



^(6>) = - - ^ . _, 2 - (59) 



7o r° 

ct± )T (0 = - rp e , e (0) / dft>ft> 3 |g(ft>,0l 2 M£(a>>"/c)- (57) 

2^(0q Jo 

In both expressions, Eqs (ISTI) and (|57l) . the spectral weighting factor |g| 2 results as a 
Lorentzian profile modulated by a damped harmonic oscillation, 

\g((0,t)\ 2 = -P(co) {1 +e- 2w -2cos[(a> - £»o>]<r w } , (58) 
7o 

where the Lorentzian spectrum is given by 

J_ % 

7T(tt>- 0^) 2 + To 

Due to the Lorentzian, the cosine in Eq. (1581) contains the argument (ft) — ft)o)? < Jot. For 
7o? < 1 we may therefore approximate Eq. (1581) by 

|g(ft),OI 2 ~-/ 5 (w)(l-^ 7o? ) 2 - (60) 
7o 

This approximation is consistent both with the fact that \g((Q, t = 0) | 2 = and with the limiting 
value for Jot ^> 1. Therefore, it may be used also for large times without introducing large 
error. 

Using the approximation (l60l) . the off-diagonal elements finally result as 

o ±iT (t) = ^A^r)(l-e-™)\ (61) 

where the spectrally averaged dissipative dipole-dipole pattern reads 

2 f°° -i 

A t( r ) = — t dcocoP(co)iJ.p(cor/c). (62) 
b ft>o Jo 

Performing analogous approximations for the diagonal matrix elements (ISTT) results in 

a ± , ± (0 = ^^A^(0)(l-e-^) 2 . (63) 
As (0) = [1 + (7o/ft>o) 2 ] 3//2 ~ 1, the latter diagonal elements become 

tf ± . ± (0^^(l-^ 7,) f- (64) 



Note, that the inner matrix elements, Eqs. (I6TI) and (l64l) . have only one time scale, y^ 1 , 
at which they increase in time. This is different from the results obtained for two atoms 
with classical fixed positions ll2TTl . There, two time scales have been derived, where the first 
coincides with y^ 1 and the second is [yo/i£ (&o r ci)] ~\ where r c \ is the classical fixed inter- 
atomic distance. 
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3.3. Averaged dissipative dipole-dipole pattern 

For the extreme far field, where the mean interatomic distance r^> r c with the critical distance 
being defined by 



(65) 



in the average (1621 the spectral width of the Lorentzian generates in the argument of the dipole- 
dipole pattern a corresponding width Ax ^> lit. Thus, the 27T -periodic oscillation of 
will be effectively averaged over to result in A« (r) ~ 0. Thus, in this case the off-diagonal 
matrix elements are 

o±, T (0 = 0, (r>r c ). (66) 

In the opposite case, when the mean inter-atomic distance satisfies the condition f < r c , 
which still permits very large distances between the atoms, the Lorentzian generates a 
sufficiently small width to provide a non- vanishing average (A*). Under these circumstances 
we may take the dipole-dipole pattern at the maximum (Oq of the Lorentzian, to obtain 



A^(r)^^(k r). (67) 

As Eq. (|67T) decays naturally when extending to the far and extreme far field, we may 
use it as a general approximation valid in the near, inductive, and far field — including the 
extreme far field. Thus, the off-diagonal matrix elements can be further approximated as 

o^{t)^ P -^ h {l-e-^f . (68) 

The average over the dissipative dipole-dipole pattern is illustrated in Fig. [TJ where the 
function (kor) is shown for the case when both atomic dipoles are perpendicular (solid 
curve) or parallel (dotted curve) to the distance vector between the atoms. As — 1 < /X* < 1, 
also the average is bound within the interval — 1 < it* < 1. But only for a strongly localized 
interatomic wavepacket that is positioned near the anti-nodes of the oscillation of and in 
the near field f ~ A, this average is non- vanishing. In all other cases, e.g. far field, strongly 
localized near a node, or weakly localized, this average vanishes. In this latter case the inner 
off-diagonal elements of the reduced electronic density matrix are <J±^(t) = for all times. 

4. Concurrence 

4.1. Concurrence of the electronic subsystem of the two atoms 
For the density matrix d36l) the concurrence takes the form 

C(t) = 2max{0,ci(0,c 2 (0}, (69) 

where 



c\(t) = \p e , g (t)\ - Ja+ i+ (t)a--(t), (70) 



Cl{t) = |0+ _(0| - JpeJt)p g Jt). (71) 
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Figure 1. Dissipative dipole-dipole pattern for atomic dipoles being perpendicular = 1, 
solid curve) and parallel (| = 0, dotted curve) to the inter-atomic distance vector. In the inset 
the dependence of the averaging over inter-atomic distances on the mean distance f and rms 
spread Ar is illustrated. 



Using the results in Eqs (1441) . (1451) . (1641) . and (1681) . these expressions can be written as 

c 1 ( £ ) = |p e!g (0)| £ -^(l- e ) 2 , 

C2(£) = ^(0)^ (M(1 _ e) 2 



(72) 



l-PfM + p eA 0) £ - 3 -p eA 0) £ 2 



E\ Pe,e(0) 



(73) 



where we defined < £ < 1 with 

e = exp(-7oO. (74) 

It can be shown that each of these functions, c\ and C2, has exactly one zero within the 
interval £ E [0, 1], £\ and £2, respectively. Moreover, these zeros are always ordered temporally 
as 



£1 > £2, 



(75) 



which means that first the concurrence vanishes at time t\ = — In £1/70 (ESD) and then revives 
at a later time t2 = — In £2/70 (ESB). Thus exactly one single ESD and one single ESB exist 
for this system. 
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This behavior differs from that obtained within the framework of a master equation where 
the atomic positions are considered as classical and fixed values [21]. In the latter case, more 
than a single ESD and ESB may occur, depending on the initial conditions (see for example 
Fig. 2 of Ref. ETI ). The difference is due to the fact that here a different regime is considered. 
Namely, we allow the atoms to move and to receive the photon recoil, whereas in Ref. |f2Tfl 
their positions are fixed (Lamb-Dicke limit). 

Returning to the concurrence functions ([721) and (f73l) . without loss of generality, the 
initial electronic state of the two atoms may be written as the reduced density matrix 

( p qs /p{\-p) \ 





P(0) 







\qy/p(\-p) l-p 
which for q = 1 reduces to the pure state 



(76) 



lv(o)) = Vp\tt) + v l -p\U)- (w) 

With this choice it can be easily seen from Eqs d72l) and (|73l) that c\ is a function of p 
and q, whereas C2 depends only on p. Furthermore, only the function C2 depends on the 
averaged dissipative dipole-dipole pattern . As a consequence the two corresponding zeros 
are functions E\ = £\{p,q) and £2 = e2(p,p.p). 

For a pure state (q = 1) the dependence of the zeros on the excited-state probability p 
and on the mean dipole-dipole pattern p^ is shown in Fig. [2] where £\ is plotted in black 
and £2 in blue (p^ = 1), green (p^ = 1/2), and red (p^ = 0). It can be observed that for a 
given value of p, the time of disentanglement, i.e. the time window between ESD and ESB, 
decreases with increasing mean p^ . 

A trivial case is obtained for p = 0, where both atoms start in their electronic ground 
states. In this case the concurrence vanishes at all times. For < p < 1, which includes 
the case of initial maximal entanglement (p = 1/2), the concurrence first decays, ending 
with an ESD, stays absent for a period of time, and then revives in the form of an ESB, see 
Fig. [3] Apart from the case p^ = (red curve) where the concurrence does not revive — 
which applies to far field, weak localization, and strong localization near a node of , the 
concurrence always revives and reaches a stationary value. Thus, there is no asymptotic decay 
to zero! 

The special case p = 1, where initially both atoms are in their electronic excited state 
and thus initially are not entangled, is of particular interest. In this case a ESB is present 
without any previous ESD nor entanglement, see Fig. HI This behavior can be explained 
when considering that entanglement between the atoms and the em field is created when one 
of the atoms emits a photon, with the identity of the emitting atom being unknown. We will 
come back to this issue in Sec. 14.21 

In the case of a mixed state, c\ and thus also the zero £1, are modified by the value of q. 
This corresponds to an obvious modification of the initial value of the concurrence together 
with a modification of the decay and thus ESD. However, the rebirth of the concurrence is 
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Figure 2. Zeros of the concurrence functions c\ (black) and ci (colored) as functions of the 
excited-state probability p. The zeros of ci are plotted for values of p.^ being: 1.0 (blue), 0.5 
(green), 0.0 (red). 



determined by the function C2 and thus is not affected, as can be seen when comparing Fig. [3] 
with Fig. [51 

4.2. Stationary concurrence due to a dark state 

For t — > °° the two functions, Eqs (1721) and (1731 reach the stationary values 

limd(0=-^, (78) 
limc 2 (0=^f^, (79) 
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Figure 3. Concurrence for an initial pure maximally entangled state (q = 1, p = 1/2) as a 
function of time for different values of jl^ : 1.0 (blue), 0.5 (green), 0.0 (red). 

so that there exists a stationary value for the concurrence, 

limC(0 = ^W (80) 

This result might at first sight appear counter-intuitive, as it states that there will be never 
any asymptotic decay to zero of the concurrence. However, when considering the electronic- 
photonic and motional dynamics it becomes clear that during the interaction of the atoms with 
the em field a so-called "dark state" is created. Starting from the two atoms being excited, 
the system may evolve into a superposition state where one photon has been emitted by either 
atom, see Fig. [6l 
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Figure 4. Concurrence for the initially separable state p = 1 (q = 1), i.e. both atoms being 
initially excited, as a function of time for different values of p,^ : 1 .0 (blue), 0.5 (green), 0.0 
(red). 

This superposition state is an eigenstate of the Hamiltonian (|4]) when neglecting the 
quantum dispersion, i.e. discarding the kinetic energy term. As a consequence, this state 
is stationary and thus a further photon emission does not occur — i.e. it is a dark state. This 
stationary state manifests entanglement between all partitions. If the emitted photon does not 
carry information on its emission center, i.e. on the identity of the emitting atom, tracing 
over the em field does not destroy the entanglement between the two atoms. Moreover, if 
the associated photon recoil does not provide information on which atom emitted the photon, 
tracing over the inter-atomic distance neither diminishes the electronic entanglement between 
atoms. In such a case, as the dark state is stationary, the concurrence reaches a stationary 



Dynamics of entanglement between two free atoms with quantized motion 



16 



c 




Figure 5. Concurrence for the initially mixed state p = 1 jl, q = 0.3, i.e. a degraded maximally 
entangled state, as a function of time for different values of p.^: 1.0 (blue), 0.5 (green), 0.0 
(red). 



value given by Eq. (l80l) . 

Contrary to this situation, when the emitted photon and/or the associated recoil carry 
information on the identity of the emitting atom, after tracing over the em field and inter- 
atomic distance, the reduced electronic state does not reveal entanglement. This situation is 
realized in various ways, all characterized by the condition p,^ = 0. For instance p,^ = if: 
(a) the atoms are weakly localized, (b) strongly localized but near a node of , or (c) in the 
far field r ^> Xq. In the latter case one may speak of two statistically independent reservoirs 
for the two atoms, where our result is in agreement with Ref Q. 
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Figure 6. Dark-state creation from state | ft) by the emission of one photon: The system 
evolves into a superposition of states | fj,) |photon) and | |photon) . 

Of course, as the kinetic-energy term in Eq. (0) has been approximated, the discussed 
darkstate is not an exact eigenstate of the full Hamiltonian d4]). As a consequence, in a 
non- approximate treatment the dispersion of the inter-atomic distance would lead to a finite 
lifetime of this darkstate. However, this lifetime would be much larger than the lifetime of the 
electronic excited state, In/jo (cf. discussion in Sec. [B- In this sense, our approximation is 
consistent for typical timescales of electronic decay. 

5. Summary and Conclusions 

In summary it has been shown that different from the case of two atoms with classical fixed 
positions, where more than one ESD and one ESB may occur ETTl . for the case of two free 
atoms with quantized relative motion and initially at rest, only one ESD and one ESB can 
occur. Furthermore, the inclusion of photon recoil permits the creation of a dark state by 
spontaneous emission. This dark state manifests itself in a stationary electronic entanglement 
between the atoms, that even may be created from an initially separable state. Whether this 
stationary entanglement appears or not, depends crucially via the photon recoil on the statistics 
of the inter- atomic distance. 

We note, that a similar behavior has been observed in the non-Markovian calculation of 
Ref. Il26ll . We believe that these phenomena can be well explained within the framework of a 
dark state, as discussed here. 
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